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ABSTRACT 

We discuss the radial dependence of the peak flux and the fluence 
of solar flare produced energetic particles under the assumption that 
they propagate diffusively in the heliosphere. 


MODEL 


There is considerable evidence that in many cases the propagation 
of solar energetic particles can be described by Parker's spherically 
symmetric transport equation which includes the effects of diffusion, 
convection, and adiabatic energy loss in the expanding solar wind 
(Parker, 1965). 
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The omnidirectional flux is then 

j - vU/4* (2) 

where v - particle speed. 

Equation (1) has generally been solved under the assumptions that 

K - K r b (3) 

r o 

where K - diffusion coefficient at 1 AU and 

o 

r - radial distance in AU. 

When the differential number density is expressed as a power law 
in kinetic energy, U - U o T’ 7 ,the explicit energy dependence of (1) can 

be eliminated. Even with these assumptions, analytic solutions to (1) 
have been found only for special cases of the radial dependence of K : 
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b - 1 (Fisk and Axford, 1968) and b - 0 (Lupton and Stone, 1973). 
Therefore (1) is usually solved numerically (Webb and Quenby, 1973; Ng 
and Gleeson, 1975; Hamilton, 1977; Zwickl and Webber, 1977). Examples 
of numerical solutions to (1) are shown in Fig. 1. 



Fig. 1. Intensity- time profiles at (a) 1 AU and (b) 5 AU from a 
numerical solution of (1) for K-1.5 x 1CT cm /s and b-0.5, 

O 

values typical for 10 MeV protons. The three curves 
correspond to diffusion only (V-0) , diffusion + convection 
(V-400 km/s , 7-I.O), and diffusion + convection + adiabatic 
deceleration (V-400 km/s, 7-3.5), (From Hamilton, 1981). 

If the effects of convection and energy loss are removed from (1) 
(by setting V-0, for example), then a pure diffusion equation 
results, and a solution has been given by Parker (1963) for any value 
of the radial index b. For the case of 3-dimensional space, 
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U(r,t) a C jr~i ) 2 " b exp C * 
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Of interest are the radial dependences of the time to maximum 
intensity, the maximum flux, and the fluence. 
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(independent of K or b) (7) 
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F - J^* j(r,t) dt (fluence) 

Fa r- (b+1) 


These analytic solutions to the pure diffusion equation are 
useful for comparison with numerical solutions to the complete 
transport equation, but generally they do not agree well with 
observations of solar energetic particles at moderate energies (< 100 
MeV) . At these energies the effects of convection and energy loss 
become important. Including these additional terms reduces t^^ and 

produces a more rapid decrease of J max and p with increasing r. 
OBSERVATIONS 

Observations of solar energetic particles beyond 1 AU have been 

made with the Pioneer 10/11 and Voyager 1/2 spacecraft. Observations 

inside 1 AU have been reported from Helios 1/2. We review here the 

two studies which used simultaneous observations at two or more radial 

distances to deduce radial propagation parameters. 

The power law index b of the radial diffusion coefficient has 

been determined for several events. Hamilton (1977) analyzed 11-67 

MeV protons in solar particle events covering the radial range 1-6 AU. 

His values for b ranged from 0.3 to 0.5 with the trend towards the 

smaller value at larger radial distances. Beeck et al. (1987) studied 

protons and heavier ions over a somewhat lower energy range (0.4-27 

MeV/nuc) for two particle events covering the radial range 0.65-1.9 

AU. Their value for b ranged from 0.5 to 0.7. In both of these 

studies, b was deduced by fitting t at two or three radial 

° max 

distances . 

To summarize these results, a value of b - 0.5 ± 0.2 covers all 
six events studied, with a trend to vary from b — 0.7 near 1 AU to 
b - 0.3 near 5 AU. Other observations at larger radial distances 



indicate a further reduction to b « 0 or somewhat less by 10 AU 
(Webber and Goeman, 1979). This variation of b with r of course 

implies that is not really a simple power law in r. Nevertheless, 

it is a useful approximation over limited radial ranges. 

From ( 6 ), we expect t a r 1,3 for b - 0.5 in the pure diffusion 

approximation, and this is close to what is observed. The radial 
dependence of t is only slightly weaker than this (Hamilton, 1977). 

.3 

Equation (7) predicts an r dependence for j . The solution to 

( 1 ), on the other hand, gives an_ r dependence for parameters 
typical for 30-67 MeV protons and r 3-3 to r ' 3 for 11-20 MeV protons 

depending on the spectral index 7 . The more rapid decrease of 

j^^with r results largely from energy loss and thus is more important 

at lower energies and for larger values of 7 . 

Hamilton(1977) also studied two events for which the maximum flux 
fell off more rapidly with radial distance (r 3 ' to r ). These 
events were observed in solar wind rarefaction regions in which the 
flux tube cross section increases more rapidly with r than « r as is 
appropriate for 3-dimensional isotropic space. Parker (1963) has 
shown this more rapid flux tube expansion results in a more rapid 
decrease of j with increasing r. 

To my knowledge, no observations of the radial dependence of the 
particle fluence have been reported. To make an estimate we may be 
guided once again by the pure diffusion approximation. Equation ( 8 ) 
would then suggest a r ‘ dependence for b - 0.5. Including effects 
of convection and energy loss causes a more rapid decrease. A 
resonable estimate is F a r 2 to r 2 ‘ 3 . 


APPLICABILITY OF RESULTS 

The model discussed above applies to a situation of spherical 
symmetry with istropic diffusion. Ng and Gleeson (1971) have shown, 
however, that the model (Eq. 1) applies anywhere within the flux tube 
connected to the flare site even though diffusion parallel to the 
interplanetary magnetic field is much more rapid than that 
perpendicular to it. 

The particle events selected for the two studies cited above were 
some of the very few observed when two or more spacecraft at different 
radial distances have simultaneously been well-connected to a flare 
site. In most events, the rapid particle intensity decrease away from 
the best connected field line causes large flux differences at 
separated spacecraft in addition to any radial dependence. Thus the 
radial dependences in J max and F cited above will rarely be observed 

in individual events except for fortuitously located spacecraft. On 
the other hand, these predictions may be useful on a statistical basis 
in extrapolating from the large data base collected at 1 AU. 

Finally, we note that there are particle events in which there is 
very little interplanetary scattering, particularly inside of 1 AU 
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(e.g. Bieber et al., 1980). These "scatter free" events are not 
described by (1). However for most of the large events in the 1-100 
HeV energy range, (1) generally appears to be a good approximation. 
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